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We study the Non-Hermitian quantum mechanics for the discrete system. This paper gives an 
exact analytic single-particle solution for an A^'-site tight-binding chain with two conjugated imagi- 
nary potentials ±47 at two end sites, which Hamiltonian has parity-time symmetry {VT symmetry) . 
Based on the Bethe ansatz results, it is found that, in single-particle subspace, this model is com- 
prised of two phases, an unbroken symmetry phase with a purely real energy spectrum in the region 
7 ^ 7c and a spontaneously-broken symmetry phase with N — 2 real and 2 imaginary eigenvalues in 
the region 7 :^ 7c. The behaviors of eigenfunctions and eigenvalues in the vicinity of 7c are investi- 
gated. It is shown that the boundary of two phases possesses the characteristics of exceptional point. 
We also construct the equivalent Hermitian Hamiltonian of the present model in the framework of 
metric-operator theory. We find out that the equivalent Hermitian Hamiltonian can be written as 
another bipartite lattice model with real long-range hoppings. 
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I. INTRODUCTION 

Since the discovery that a non-Hermitian Haniiho- 
nian having simultaneous parity-time (VT) symmetry 
has an entirely real quantum-mechanical energy spec- 
trum fl|, there has been an intense effort to establish a 
T^T-symmetric quantum theory as a coinplex extension 
of the conventional quantum mechanics [3, H, |3, H, H, 0] ■ 
Although there have been no experiments to show clearly 
and definitively that quantum systems defined by non- 
Hermitian T'T-symmetric Hamiltonians do exist in na- 
ture, many models have been proposed to verify theorems 
and perform numerical and asymptotic analysis (for a re- 
cent review, see [s] and references therein). The reality 
of the spectra is responsible to the VT symmetry: If all 
the eigenstates of the Hamiltonian are also eigenstates of 
VT, then all the eigenvalues are strictly real and the sym- 
metry is said to be unbroken. Otherwise, the symmetry 
is said to be spontaneously broken. In practice, imag- 
inary potential usually appears in a system to describe 
physical processes phenomenologically, which have been 
investigated under the Non-Hermitian quantum mechan- 
ics framework 

SEIiilMIilBEIiiliilill- how- 

ever it is not clear whether non-Hermitian Hamiltonians 
describe real physics or are just unrealistic mathemat- 
ical objects. It is known, for a diagonalizable Hamil- 
tonian, that the presence of PT-symmetry implies the 
pseudo-Hermiticity of the Hamiltonian [i9|, i.e., VT- 
symmetry is a special case of pseudo-Hermiticity. A di- 
rect way of extracting the physical meaning of a pseudo- 
Hermitian Hamiltonian having a real spectrum is to seek 
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for its Hermitian counterparts [l^ [3 • The metric- 
operator theory outlined in 0] provides a mapping of such 
a pseudo-Hermitian Hamiltonian to an equivalent Hermi- 
tian Hamiltonian. The construction of the latter is usu- 
ally quite complicated if the Hilbcrt space of the systems 
is infinite-dimensional. It is generally more tractable for 
lattice systems with a finite-dimensional single-particle 
Hilbert space. Under such circumstances the study of 
simple models which are exactly solvable and slight mod- 
ified version of a practical model, and yet are tractable, 
is of particular importance to examine the impact of an 
imaginary potential on the feature of eigenvalues and 
eigenfunctions of a lattice system. 

In this paper, we will focus on the non-Hermitian VT- 
symmetric quantum theory for a discrete system [20I I2H , 
the original form of which is exploited to describe the 
solid-state system in condensed matter physics or cou- 
pled quantum devices since the advent of quantum in- 
formation theory [l^j- As an illustration, a simple iV- 
site tight-binding chain with the uniform nearest neigh- 
bor (NN) hopping integral is concerned. Such a model 
is used to describe the Bloch electronic system in con- 
densed matter physics and now the qubit array relevant 
to quantum information applications. The correspond- 
ing non-Hermitian T^T-symmetric version is constructed 
by adding two conjugated imaginary potentials ±17 at 
the end sites. The objective of this paper aims at the 
exact solutions of such a model so as to confirm the non- 
Hermitian T^T-symmetric quantum theory for the dis- 
crete system. Based on the Bethe ansatz results, it is 
found that, in single-particle subspace, this model ex- 
hibits two phases, an unbroken symmetry phase with a 
purely real energy spectrum when the potentials are in 
the region 7 7c and a spontaneously-broken symmetry 
phase with N ~2 real and 2 imaginary eigenvalues when 
the potentials are in the region 7 )^ 7c Based on the exact 
solutions, the behaviors of eigenfunctions and eigenvalues 
in the vicinity of 7c are investigated. It is shown that the 
boundary of two phases possesses the characteristics of 
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exceptional point: there are two eigenstates coalescing at 
7c as well as with square-root type level repulsion in the 
vicinity of it. We also construct the equivalent Hermi- 
tian Hamiltonian of the present model in the framework 
of metric-operator theory. We find out that the equiv- 
alent Hermitian Hamiltonian can be written as another 
bipartite lattice model with real long-range hoppings. 

This paper is organized as follows, in Sec. II, the model 
is presented. In Sec. Ill, the Bethe ansatz solutions are 
given. In Sec. IV, we investigate the characteristics of the 
critical point 7^. Sec. V is devoted to the construction 
of the equivalent Hermitian Hamiltonian. Sec. VI is the 
summary and discussion. 



II. PT-SYMMETRIC UNIFORM 
TIGHT-BINDING CHAIN 

We consider a simplest discrete system described by a 
non-Hermitian Hamiltonian H having the VT symmetry. 
It is a tight-binding chain with uniform nearest neighbor 
hopping integral and two additional conjugated imagi- 
nary on-site potentials on the two end sites, which can 
be written as follows: 



H 



N-l 



a\ai+i +'ii.c.\ +i^a\ai- i^a\aN, (1) 



where a] is the creation operator of the boson (or 
fermion) at lih site, the tunnelling strength and poten- 
tial are denoted by J and ±^7. V and T represent 
the space-reflection operator, or parity operator and the 
time-reversal operator respectively. The effects of V and 
T on a discrete system are 



TiT = -i, Va]V 



a 



t 

N+l-l- 



(2) 



Obviously, the Hamiltonian Eq. ([T]) has VT symme- 
try, i.e., H'^'^ = VTHVT = H. According to the 
PT-symmetric quantum mechanics Q , H can be further 
classified to be either unbroken VT symmetry or bro- 
ken VT symmetry, which depends on the symmetry of 
the eigenstates jV'/c) in different regions of 7. The time- 
independent Schrodinger equation is 



(3) 



with corresponding eigenvalue e^. The system is unbro- 
ken VT symmetry if all the cigenfunctions have VT sym- 
metry 



VT\^k) = \i>k) 



(4) 



and all the corresponding eigenvalues are real simulta- 
neously. This classification depends on the value of the 
parameter 7. Beyond the unbroken VT symmetric region 
the system is broken VT symmetry, where Eq. Q does 
not hold for all the cigenfunctions and the eigenvalues of 



broken VT symmetric cigenfunctions are imaginary. One 
of the aims of this paper is to provide the complete exact 
cigenfunctions, eigenvalues, and the boundary between 
unbroken and broken VT symmetric regions. 

Acting the VT operation on an arbitrary single- 
particle wave function 

\^) =Y,hia\\Q) , (5) 

where |0) denotes the vacuum state, we have 

VT 1^) - e^')* «lv+i-/ 10) - E i^N+i-iT A |0) . 

(6) 

Then in the rest of this paper, we can simply use 

VThi = (7) 

to present the VT operation on the single-particle wave- 
function. 

III. BETHE ANSATZ SOLUTIONS 

In this paper, we only focus on the single-particle case. 
We denote the single-particle eigenfunction in the form 



IV'fe) = fia\ |0) 



(8) 



The Bethe ansatz wavefunction can be expressed in the 
form 



fi = A{k)e''''+B{k)e 



ikl 



(9) 



where the quasi- momentum k and amplitudes A(k), 
B (fc) can be determined from the Eq. ([3]) and the proper 
definition of the inner product according to the VT- 
symmetric quantum theory Q. The solutions are pre- 
sented explicitly in the following. 



A. Unbroken "PT-symmetric region 

In the unbroken VT symmetric region 7 ^ 7c where 



7e= + l 

{ J, N ^ 2n 

the explicit CVT normalized wavefunctions are 



^ik{l~No) 



T] (fc) e 



-ik(l+No) 



l + \?]{k)\^ sin {Nk) I sink- 2Nr] (fc) e-^^^^+i) 

(11) 
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where Nq = {N + 1) /2 is the center of the chain and the 
coefficient 



76 — iJ 



(12) 



■ye — iJ 

The quasi-momentuni k satisfies the equation 

7^ sin [k{N - 1)] + J2 sin [fc(iV + 1)] = (13) 

which has N real solutions. All the corresponding eigen- 
values are real 

Efc = -2Jcosfc, (14) 
with the quasi-momentum k being more explicit form 



k = — ,nk e [l,iV] 



and 



9fc = tan ^ 



(7^ - J') 
(7^ + J2) 



tan k 



(15) 



(16) 



The reality of the spectrum is a consequence of VT in- 
variance. At 7 = 0, we have 77 (fc) = 1, the eigenfunctions 
reduce to the form 



fth = 0) = i-^r^ 



N+1 



sm{kl), (17) 



k - -^^,nk e [1,N] 



N + 1' 

which is the well-known solution of a Hermitian tight- 
binding chain. 

B. Broken PT-symmetric region 

In the region 7 7c, the VT symmetry of the Hamilto- 
nian is spontaneously broken; even though VT commutes 
with H, the eigenfunctions of H are not all simultane- 
ously eigenfunctions of VT. In this region, it can be 
shown that there are N — 2 real k for the equation Eq. 

which corresponds to the eigenfunctions Eq. (fTTI) 
and real eigenvalues Eq. dH]). The rest two eigenfunc- 
tions correspond complex quasi-momenta, which are in 
the form k = 7r/2 ± i \k\ with k satisfying the equation 

72 sinh [k{N - 1)] = J2 sinh [k{N + 1)] , (odd N) ; 
7^ cosh [k{N - 1)] = J2 cosh [k{N + 1)] , (even N) . 

(18) 

The broken T^T-symmetric eigen functions can be written 
as 



/2±iK 



oc e 



±K.No 



J 



J + je 



±K 



with the imaginary eigenvalues 

Ek = ±i2 Jsinh n. 



(19) 



(20) 



At this stage, the wave functions Eq. (|T9|) are not nor- 
malized since the standard CVT normalization procedure 
is invalid for such broken PT-symmetric wave functions. 



C. CVT formalism 

In this subsection, we will elucidate the CVT formal- 
ism based on the solutions of the present model. It can be 
seen that the solutions have different features compared 
to that of a Hermitian Hamiltonian. The most intuitive 
novelty of the eigenfunctions Eq. (|lip with real eigen- 
values are not standing waves due to |?7(fe)| being not 
unitary for nonzero 7. According to the PT-symmetric 
quantum theory 7], the correct inner product is deter- 
mined by the Hamiltonian itself, so it is necessary for 
the discrete system Eq. ([T]) to establish a self-consistent 
formalism based on the obtained exact solutions. 

A straightforward calculation shows that the eigen- 
functions Eq. pT|) obey the T'T-symmetry 

VTfl - fl (21) 

According to the PT-symmetry quantum theory 0] , one 
can define the coordinate-space representations of the 
parity operator V and C operator for a discrete system 
as 



V{m,l) 

C{m,l) = Y^f^f^ 



I 



(22) 
(23) 



which lead to the CVT orthogonal and normalization re- 
lations 



(24) 



l,rn 



This guarantees the positive-definite CVT inner product 
of two arbitrary states and that the time evolution is 
unitary. One can verify that the operator C satisfies 



= 1, [C,VT] = 0, [C,H] = 0. 



(25) 



In the next section, the above formalism will be uti- 
lized to establish the canonical basis to construct the 
equivalent Hermitian Hamiltonian in the unbroken VT- 
symmetric region. We will see that choosing CVT nor- 
malized eigenstates of H leads the Hermitian equivalence 
matrix to be more symmetrical. 



IV. EXCEPTIONAL POINTS 

In this section, we investigate the critical behavior of 
the system as 7 in the vicinity of 7c. From the obtained 
solutions of the model, the critical point 7c has the char- 
acteristics of exceptional point [23, 24, 25, 26]. We will 
investigate the feature of eigenvalues and eigenfunctions 
around the exceptional point in detail. 

In unbroken T'T-symmetric region, as 7 approaches to 
7c, the solutions of Eq. change abruptly. Actually, 
as the function 

jr(fc) = 72 sin [k{N - 1)] - J2 sin [k{N + 1)] (26) 
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is an even (odd) funct 
(odd) N, two real energy 
disappear when 7 passes 
two regions. Meanwhile, 
pear. Thus the critical 
be characterized by the 
7 = 7c — 0+, we denote 



fl/2+5 and 4 



ion about k = Tr/2 for even 
levels, which are closest to zero, 
through the boundary 7c of the 
two imaginary energy levels ap- 
behavior of the eigenstates can 
two components problem. For 
the referring eigenfunctions as 



/2-S 



with 




N is even 



iVisodd 



r 



J2' 



(27) 
(28) 



It follows that such two eigenfunctions approach to a 
same function and their VT norms tend to zero when 
VT symmetry is broken. 

In the broken PT-symmetric region, for 7 = 7c + 0"*", 
such two eigenfunctions are replaced by Eq. (jl9p . while 
the corresponding eigenvalues turn to imaginary values 
Eq. (PO)) . Other N — 2 eigenfunctions in the form of Eq. 
(fTTjl with real eigenvalues still satisfy the CVT orthog- 
onal and normalized relations Eq. ([24]) . However, the 
two eigenfunctions Eq. (|19p can no longer be normalized 
via the above CVT inner product since they have the 
following features 



VT 



/7r/2+i 



a li/2- 



and 



J TT 



f 

J 71 



/2±i 



= 0. 



(29) 



(30) 



On the other hand, the corresponding energy levels 
experience a switch from real to complex values as well 
as a coalescence at the exceptional point. In the unbroken 
symmetric side, from (|27l) we have 

e^/2±s ^ ±2 J sin 6 (31) 
~ ±2JS. 

In the broken symmetric side, from (|18p we have 

N is even 
TV is odd ' 

which correspond to the eigenvalues 

£^/2±iK - ±2iJsinhK 
~ ±2ijK. 




(32) 



(33) 



Alternatively, taking I7 — 7c| as the variable, we can see 
that the two concerned eigen states satisfy 



£Tr/2±S — Ini(e7r/2±iK) 



±2jJ^^ TV is even 
±2jJ%^, TV is odd 



(34) 



near the critical point, which reveals the symmetry of 
the critical behavior. In Fig. ^ we plot the real and 
imaginary parts of two repelling levels as functions of 
"f — jc, which are obtained from the approximate ana- 
lytical results Eqs. ([3T|) . p3|) and numerical simulations 
for finite systems. The analytical eigenvalue expressions 
(I3ip and p3p are good approximations to the numeri- 
cally exact results. It shows that the exceptional points 
are always associated with a level repulsion in the vicinity 
of them. The square-root type functions for the energy 
reveal this characteristics. Previous study [U,!!^, shows 
that the exceptional point is often related to the emer- 
gence of chaotic behavior. However, the quantum chaos is 
not found in the present model. This may be due to that 
the coalescence and repulsion of levels in this model only 
relate to two eigenstates rather than multi-level. Thus 
there is no occurrence of quantum chaos. 



V. EQUIVALENT HERMITIAN HAMILTONIAN 

A natural question to ask is whether such a model Eq. 
(H]) describes real physics or is just an unrealistic mathe- 
matical product. While there is as yet no answer to this 
question, we can gain some insight regarding its equiv- 
alent Hermitian counterpart 0, This section aims 
at seeking the equivalent Hermitian Hamiltonian of the 
present model in the framework of metric-operator the- 
ory p. 

We start with the eigenstates of the Hamiltonian W = 
H* = H{-j). From the explicit CVT normahzed 
wavefunctions of within the unbroken PT-symmetric 
region are 



Jk(l-No) 



C(fc)e- 



-ik{l+Na) 



1 + IC (fc) 1^ sin (Nk) / sin fc - 2NC (k) e-»'=(^+i) 

(35) 



where the coefficient is 



C(fc) 



ik 



iJ 



(36) 



and the corresponding eigenvalues are (|14p . the same as 
that of Within the unbroken T'T-symmetric region, 
we have 



VT\gk)^\gk) 



(37) 



and two sets of eigenfunctions {gfc} and {/^} form a 
biorthonormal system, i.e.. 



(38) 



According to the metric-operator theory 0], one can 
construct a positive-definite operator 
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-6 -4-2 2 
y-y (unitsofj) ^^q" 




2 4 6 

y-y (unitsofj) 



o) -0.005 




2 4 6 8 10 

y-y^(unitsof J) ^ ^g-i 



v+ = \9k) {gk\ 



(39) 



From the Appendix, it can be shown that this operator 
satisfies 



VTr]+VT = r/+. 



(40) 
(41) 
(42) 



In the spatial coordinate space spanned by basis |m) 
alyi |0), the matrix representation of 77+ has the form 



(43) 



More exphcitly, it can be shown that the matrix has the 
foUowing properties 



(m| 77+ \n) = {N +1- n\ 77+ |iV + 1 - m) . (44) 



and 



(mh+|n)^(~l)"+"(m|,y+|n)*. (45) 
In addition, defining the canonical transformation R, 



we have 



(46) 



(47) 



These features allow characterizing the equivalent physi- 
cal system of -ff . In the eigenspace of 77+ , we can rewrite 
ri+ as 



E 



En en) (Cn 



(48) 



where |e„) is the eigenvector of operator 77-|_ with the 
eigenvalue £«, i-e.. 



FIG. 1; (Color online) Real and imaginary parts of two re- 
pelling levels as functions of 7 — 7c. The plots are obtained 
from the approximate analytical results Eqs. (|33|l and 

numerical simulations for the systems with A'' — 19, 20, 199, 
and 200. It shows that the analytical eigenvalue expressions 
are good approximation to the numerically exact results, es- 
pecially for small N system. 



?7+ |e„) = e„ |e„) . (49) 

The eigenvalues e„ are all real due to the Hermiticity of 
77+ and the complete set {|en)} is referred as a canonical 
metric basis. 

Accordingly, the equivalent Hermitian Hamiltonian Ti. 
can be obtained by a unitary transformation and can be 
expressed as 



Ti- — ^ y —H„in \em) (£11 



(50) 
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where Hmn = (cml H \ e„) is the matrix representation of 
Hamiltonian H under this canonical metric basis {|e„)}. 
The Hermitian equivalence matrix Ti can be achieved as 



Tin 



Then real symmetric matrix Ti can be regarded as a 
single-particle matrix representation of a tight-binding 
model on a bipartite lattice N ~ Na + Nb with 
\Na - iVsl = 1 (0) for odd (even) N. Two sublattices are 
connected by the long-range hopping terms with strength 
(51) ^ij ■ '^^^ corresponding tight-binding Hamiltonian can be 



2J 



where eL, = (z I 

^m) is the component of state Itm) i^ the 
spatial coordinate basis 

Before we explore the features of the equivalent Her- 
mitian Hamiltonian it is worthy to point that the or- 
dinal tight-binding chain (jT]) is a bipartite lattice which 
can be separated into A and B sublattices with sites Na 
and Nb^ respectively. We have N = Na + Nb with 
\Na - A''b| = 1 (0) for odd (even) N and the two sub- 
lattices are connected by hopping terms with strength J 
(see Fig. [5] (a) and (c)). Now we turn to the properties 
of TL. It can be shown, from the Appendix of this paper 
[27| . that under a proper choice of the transformation of 
{|e„)}, the matrix representation of Ti. can be in the form 
of 



A 
A^ 



(52) 



where is the transposed matrix of a N a x Nb matrix 
A, and A satisfies 



A, 



^NA+\-j,NB+\-t-. ^ is even 
^TVA+i-i.A^B+i-ji is odd 



(53) 



J 



expressed as 



1=1 j 



I AB 



(jI+H.c), 



(54) 



where \l) a \l) -g denote single-particle states referring 
to Zth site in sublattices Na and Nb, respectively. Note 
that both the original model H and its Hermitian coun- 
terpart Ti are all bipartite. The former contains NN cou- 
plings, while the latter contains the long-range couplings. 
This fact agrees with the observations from other exam- 
ple Hamiltonians, that in general the Hermitian counter- 
part of a pseudo- Hermitian Hamiltonian obtained by the 
metric-operator theory is nonlocal operator 0. 

In order to exemplify the above analysis, we consider 
the small systems with N = 7 and 8. Fig. [5] shows the 
schematics of configurations for (a,c) pseudo-Hermitian 
system H on 7, 8-site lattices and (b,d) their Hermitian 
counterparts Ti. The corresponding hopping constants 
Ay are computed for 7 — 0.00, 0.50, and 0.99, which are 
listed in Table 1 (a) and (b). It indicates that all the 
constants vary within a narrow range without changing 
signs as 7 covers the whole unbroken symmetric region 
and obey the relation of Equation ([55]) . 



Table 1 (a) 



7 


All, A34 A12, A33 


A13, A32 


Al4, A31 


A21, A24 


A22, A23 


0.00 


0.6242 1.0068 


-0.2997 


0.0830 


0.2071 


-1.2071 


0.50 


0.5703 0.9731 


-0.3089 


0.0883 


0.2039 


-1.2075 


0.99 


0.3355 0.8949 


-0.3280 


0.0774 


0.1468 


-1.2089 



Tabic 1 (b) 



7 Aii,A44 Ai2,A34 Ai3,A24 Ai4 A2l,A43 A22,A33 A23 A3i,A42 A 



32 



A4 



0.00 0.5627 -0.9300 -0.2994 -0.1199 0.1954 1.1615 -1.2411 0.0914 0.3333 0.0277 
0.50 0.5153 -0.8918 -0.3057 -0.1304 0.1909 1.1527 -1.2469 0.0972 0.3461 0.0310 
0.99 0.1766 -0.9005 -0.3157 -0.1458 0.1005 1.0333 -1.3522 0.0627 0.3505 0.0143 



Table 1. The coupling constant distributions Xij for systems with (a) N — 7 and (b) 8, obtained from numerical 
simulation for 7 = 0.00, 0.50, and 0.99. 



VI. CONCLUSION AND DISCUSSION 



In conclusion, we have studied the Non-Hermitian 
quantum mechanics for the discrete system. The ex- 
act analytic single-particle solution for a tight-binding 
chain with two end imaginary potentials is obtained, 
which substantiates the formalism of the quantum 



theory for discrete system: Based on the Bethe ansatz 
results, it is found that, in single-particle subspace, this 
model exhibits two phases, an unbroken symmetry phase 
with a purely real energy spectrum in the region and 
a spontaneously-broken symmetry phase with a pair of 
imaginary eigenvalues in the region. 

In addition, we have indicated boundary of two phases 
possesses the characteristics of exceptional point. We 
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FIG. 2: (Color online) Schematic illustration of configurations 
for (a,c) pseudo-hermitian system H on 7, 8-site lattices and 
(b,d) their hermitian counterparts h. Both the original system 
and its hermitian counterpart are all bipartite graphs. 

also construct the equivalent Hermitian Hamiltonian of 
the present model in the framework of metric-operator 
theory. We find out that the equivalent Hermitian Hamil- 
tonian can be written as another bipartite lattice model 
with real long-range hoppings. 

It is worthwhile to note that the Non-Hermitian 
Hamiltonian ([T]) cannot simply be written as the form 
^j.efca[.afc, where = 'l2i fi'^l creation operator in 
k space, since the transition matrix is no longer unitary. 
However, the equivalent Hermitian Hamiltonian Ti. can 
be in such form. It might be interesting to explore the 
multiparticle sector of the system using the equivalent 
Hermitian Hamiltonian. 

We acknowledge the support of the CNSF (Grants No. 
10874091 and No. 2006CB921205). 



PROPERTIES OF MATRICES 

77+ AND n 



VII. APPENDIX: 



According to the defination of the metric operator 
and the properties of eigenfunctions, equations ([57]) and 
we have 



V A; / k' 



(55) 



and 



VTr,+VT = VT}_^\gk){gk\VT 

k 

= E 1^'=^ ^5fel = 



(56) 



On the other hand, equation allows the real matrix 
representation of the metric operator by the transforma- 
tion on the spatial coordinate basis 

mod[;,2] 
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For the sake of simplicity, hereafter 77+ denotes a real 
matrix. Then from equations ([1T|) and (|^^ indicate 



-1 



(58) 
(59) 
(60) 



i.e., r?+ is a real unitary and bisymmetric matrix. Then 
the eigenfunctions |e„) of 77+ can always be written as 
real functions and obey 



Vie,. 



(61) 



These properties is helpful for characterizing the equiva- 
lent physical system of H. 

Now we focus on the case of even N. A straight- 
forward calculation shows that eigenfunctions |e„) and 
V |e„) have different parities. Then we can reorder the 
canonical metric basis as 

|e„) =i?|ew+i_„), ne [l,7V/2], (62) 

which satisfy 

V\en) = |e„), ne[l,N/2], (63) 
V\en) = -|e„), ne [N/2 + 1, N/2]. 

Under such canonical metric basis, we have 




(64) 



emWTHVTler. 



)3(m)+)?(n) 



-1) 



where 



i?(to) 



i5(m)+i9(n) o,* 



0, TO < N/2 

1, TO > N/2 



(65) 



is the Heaviside step function. On the other hand, since 
under the new definition of basis |^) Hmn is always imag- 
inary, all the elements Timn with d{m) -I- d{n) = 0,or 2 
must vanish. Because |e„) can always be written as real, 

one can simply take |e„) —>■ (\/— T)''^"'' jen) to obtain the 
final form of the equivalent Hermitian Hamiltonian 



n 



A 





(66) 



10. 



(57) 



where A is real matrix and A^ is its transposed matrix. 
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